Scaling of Traction Forces with Size of Cohesive Cell Colonies 
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To understand how the mechanical properties of tissues emerge from interactions of multiple cells, 
we measure traction stresses of cohesive colonies of 1-27 cells adherent to soft substrates. We find 
that traction stresses are generally localized at the periphery of the colony and the total traction 
force scales with the colony radius. For large colony sizes, the scaling appears to approach linear, 
suggesting the emergence of an apparent surface tension of order 10 -3 N/m. A simple model of the 
cell colony as a contractile elastic medium coupled to the substrate captures the spatial distribution 
of traction forces and the scaling of traction forces with the colony size. 



PACS numbers: 87.17.Rt, 87.19.R-, 68.03.Cd 

Tissues have well defined mechanical properties such 
as elastic modulus pQ. They can also have properties 
unique to active systems, such as the homeostatic pres- 
sure recently proposed theoretically as a factor in tumor 
growth [2]. While the mechanical behavior of individual 
cells has been a focus of inquiry for more than a decade 
[3HS] , the collective mechanics of groups of cells has only 
recently become a topic of investigation [7HT5]: it is un- 
known how collective properties of tissues emerge from 
interactions of many cells. 

In this Letter, we describe measurements of traction 
forces in colonies of cohesive epithelial cells adherent to 
soft substrates. We find that the spatial distribution and 
magnitude of traction forces are more strongly influenced 
by the physical size of the colony than by the number of 
cells. For large colonies, the total traction force, T , that 
the cell colony exerts on the substrate appears to scale 
as the equivalent radius, R, of the colony. This scaling 
suggests the emergence of a scale-free material property 
of the adherent tissue, an apparent surface tension of or- 
der 10 _3 N/m. A simple physical model of adherent cell 
colonies as contractile elastic media captures this behav- 
ior. 

To measure traction stresses that cells exert on their 
substrate, we used traction force microscopy (TFM) [T6] . 
Our TFM setup consisted of a film of highly elastic sili- 
cone gel (Dow Corning Toray, CY52-276A/B) with thick- 
ness h s = 27 /im on a rigid glass coverslip (Fig. [T]A). Us- 
ing bulk rheology, we estimated the Young's modulus of 
the gel to be 3 kPa. To quantify the gel deformation dur- 
ing our experiments, our substrates contained two dilute 
layers of fluorescent beads (radius 100 nm, Invitrogen): 



one layer between the glass and gel and a second at height 
z Q = 24 /im above the coverslip [17] • To image the flu- 
orescent beads, we used a spinning-disk confocal micro- 
scope (Andor Revolution, mounted on a Nikon Ti Eclipse 
inverted microscope with a 40 x NA 1.3 objective). Af- 
ter determining bead positions using centroid analysis in 
Matlab [18 , we calculated the substrate displacement, 
uf(r, 2 Q ), across ^s stressed (with cells) and unstressed 
(with cells removed) states. In Fourier space, the in-plane 
displacement field is related to the traction stresses at the 
surface of the substrate via linear elasticity, crf z (k, h s ) = 
Qij(k, z , h s )Uj(\s., z ), where k represents the in-plane 
wave vector. Here, cr f z (k, h s ) and ^|(k, z Q ) are the 
Fourier transforms of the in-plane traction stress on the 
top surface and the displacements just below the sur- 
face, respectively. The tensor, Q, depends on the thick- 
ness and modulus of the substrate, the location of the 
beads, and the wave vector [171 [19]. We calculated the 
strain energy density, w(r) = ^<jf z (r,h s )uf(r,h s ), which 
represents the work per unit area performed by the cell 
colony to deform the elastic substrate [20] . The displace- 
ment of the surface was determined using w|(k, h 8 ) = 
Qi^Qt, h 8 , h s )Q jk (K Zo, h s )u s k (k, z Q ). 

Primary mouse keratinocytes were isolated and cul- 
tured as described in [21 . We plated keratinocytes on 
fibronectin-coated TFM substrates. After the cells pro- 
liferated to the desired colony sizes over 6-72 hr, we raised 
the concentration of CaCl2 in the growth medium from 
0.05 mM to 1.5 mM. After 18-24 hr in the high-calcium 
medium, cadherin-based adhesions formed between adja- 
cent keratinocytes, which organized themselves into cohe- 
sive single-layer cell colonies [22, 23 . After imaging the 
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beads in their stressed positions, we removed the cells 
by applying proteinase K and imaged the beads in their 
unstressed positions. 

Stress fields and strain energy densities for represen- 
tative colonies of one, two, and twelve keratinocytes are 
shown in Fig. [I] Traction stresses generically point in- 
ward, indicating that the colonies are adherent and con- 
tractile. Regions of high strain energy appear to be lo- 
calized primarily at the periphery of the single- or multi- 
cell colony. For single cells, these findings are consistent 
with myriad previous reports on the mechanics of iso- 
lated, adherent cells [24-27]. Recent reports have also 
observed localization of high stress at the periphery of 
small cell colonies on micropatterned substrates [28] and 
at edges of cell monolayers [7J [8] [13] . To visualize cell- 
cell and cell-matrix adhesions, we immunostained multi- 
cell colonies for E-cadherin and zyxin. Additionally, we 
stained the actin cytoskeleton with phalloidin. Actin 
stress fibers were concentrated at colony peripheries and 
usually terminated with focal adhesions, as indicated by 
the presence of zyxin at the fibers' endpoints. In contrast, 
E-cadherin was localized at cell-cell junctions, typically 
alongside small actin fibers. Despite differences in the ar- 
chitecture of the relevant proteins, the stresses and strain 
energy distributions are remarkably similar in the single- 
cell and multi-cell colonies. 

To explore these trends, we measured traction stresses 
of 45 cohesive colonies of 1-27 cells. For each colony, we 
defined an equivalent radius, R, as the radius of a disk 
with the same area. The equivalent radii ranged from 
20 to 200 jam. We calculated the average strain energy 
density as a function of distance, A, from the colony 
edge (Fig. [5] inset). Figure [2] shows the normalized strain 
energy profiles, w(A)/w(0), of all 45 colonies. Usually, 
the strain energy density was largest near the colony edge 
(A = 0). Because of the finite spatial resolution of our 
implementation of TFM, we measured some strain energy 
outside colony boundaries (A < 0). 

Next, we examined how global mechanical activity 
changes with the cell number and geometrical size of the 
colony. As in previous studies, we calculated the "total 
traction force" [29l 130] . 

?=J W^) 2+ fe) 2 ' (i) 

exerted by the cell colony onto the substrate. This quan- 
tity is meaningful when stresses have radial symmetry 
and are localized at the colony edge, which is the case 
for the majority of colonies in this study. We observed a 
strong positive correlation between equivalent radius and 
total force over the range of colonies examined (Fig. [3]). 
Similar trends have been seen for isolated cells over a 
smaller dynamic range of sizes [29 -32 . We see no system- 
atic differences in T for colonies of the same size having 
different numbers of cells, suggesting that cohesive cells 
cooperate to create a mechanically coherent unit. 
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FIG. 1: (color online) Traction stresses and strain energies 
for colonies of cohesive keratinoctyes. (A) Schematic of ex- 
perimental setup (not to scale) with a cell colony adherent to 
an elastic substrate embedded with two dilute layers of flu- 
orescent beads. (B, D, F) Distribution of traction stresses, 
<Ji Z , and (C, E, G) strain energy, w, for a representative sin- 
gle cell, pair of cells, and colony of 12 cells. Traction stress 
distribution is overlaid on a DIC image (B) or images of im- 
munostained cells (D, F). Solid lines in (B-C, E, G) mark cell 
boundaries. For clarity, only one-quarter of the calculated 
stresses are shown in (B, D) and one-sixteenth of the stresses 
in (F). Scale bars represent 50 /mi. 

The data in Fig. [3] while scattered, show clear mono- 
tonic growth of the mechanical output of cell colonies 
with their geometrical size, independent of the number 
of cells. For smaller colonies (R < 60 /im), the increase of 
total force is superlinear. As the cell colonies get larger, 
the scaling exponent appears to approach unity. We hy- 
pothesize that the transition to an apparently consistent 
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FIG. 2: (color online) Spatial distribution of strain energy 
for colonies of different size. Each solid curve represents a 
colony's average strain energy density as a function of distance 
from the edge of the colony, A. For clarity, the profiles are 
spaced vertically according to the size of the colony. Each 
profile terminates at the point where the inward erosion of 
the outer edge covers the entire area of the colony, at A w R. 
The erosion proceeds in discrete steps of size 5, as illustrated 
in the inset. 



exponent for the large colonies reflects the emergence of 
a scale- free material property of an adherent tissue, de- 
fined by the ratio J 7 /(2ttR) = (8 ± 2) x 10" 4 N/m, with 
dimensions of surface tension. 

Just as intermolecular forces yield the condensation of 
molecules into a dense phase, cohesive interactions be- 
tween cells, mediated by cadherins, cause them to form 
dense colonies [33 , 34 . For large ensembles of molecules, 
molecular cohesion creates a free energy penalty per unit 
area, known as surface tension, for creating an interface 
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FIG. 3: (color online) Mechanical output of keratinocyte 
colonies versus geometrical size. Total force transmitted to 
the substrate by the cell colonies, defined in Eq. is plotted 
as a function of the equivalent radius of the colonies. The 
dashed line represents the scaling expected for surface tension, 
T ~ R. The solid line shows a fit of the data to the minimal 
contractility model in Eq. 



between two phases. It is tempting to think of the adher- 
ent colonies studied here as aggregates of cohesive cells 
that have wet the surface [35]. Indeed, when matter wets 
a surface, the traction stresses are localized at the con- 
tact line [36], as we found in our cell colonies (Figs. [I] 
and [2). 

Effective surface tension of cell agglomerates has been 
invoked to explain cell sorting and embryogenesis [37] . 
Previous measurements of non-adherent aggregates of co- 
hesive cells reported effective surface tensions between 2 
and 20mN/m [38ti4Q] . However, the origins of the effec- 
tive surface tension of cohesive cells are distinct from con- 
ventional surface tension. Recently, it was suggested that 
the surface tension is not only determined by contribu- 
tions from cell-cell adhesions but also the contraction of 
act o- myosin networks [4TJ 02] . It is important to distin- 
guish the effective surface tension due to active processes 
from the familiar surface tension defined in thermody- 
namic equilibrium. 

To elucidate the origins of an effective surface ten- 
sion in these experiments, we consider a minimal model 
proposed recently to describe cell-substrate interactions 
[43] [44] . We describe a cohesive colony as an active elas- 
tic disk of thickness h and radius R (Fig. [l}^). The me- 
chanical properties of the cell colony are assumed to be 
homogeneous and isotropic with Young's modulus E and 
Poisson's ratio v. Acto- myosin contractility is modeled 
as a contribution to the local pressure, linearly propor- 
tional to the chemical potential difference, A/i, between 
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ATP and its hydrolysis products [45]. In our model, 
the strength of cell-cell adhesions is implicitly contained 
in the material parameters of the colony, E and v. The 
constitutive equations for the stress tensor, cr^-, of the 
colony are then given by 



E 



2(1 + 



2v 



\-2v 



diUj + djUi + 5ij(A{i, 
(2) 

where u is the displacement field of the cell colony and 
( > a material parameter that controls the strength 
of the active pressure, £A/i. Mechanical equilibrium re- 
quires that djdij = 0. 

We use cylindrical coordinates and assume in-plane 
rotational symmetry. The top surface is stress-free, 
a rz \ z= h-\-h s = 0, and we employ a simplified coupling of 
the colony to the substrate. Ignoring all nonlocal effects 
arising from the substrate elasticity, cr r z\z=h 3 = Yu r (z = 
h s ) « Yu r . Here, u r is the radial component of the dis- 
placement field, the bar denotes z-averaged quantities, 
and the rigidity parameter, Y, describes the coupling of 
the contractile elements of the colony to the substrate. 
The local proportionality of stress and displacement is ac- 
curate only when the substrate thickness is much smaller 
than the characteristic length scale of the stress distri- 
bution or when the cells are on substrates of soft posts 

With these assumptions, the equation of force-balance 
simplifies to 



[d r (ra rr ) 



9] jr = Yu r /h. 



(3) 



Combining Eqs. Q and ([3|, we find the governing equa- 
tion for the radial displacement, u r : 



r 2 d 2 r u r + rd r u r - (l + r 2 /Q u r = 0, 



(4) 



where the penetration length, £ p , describing the localiza- 
tion of stresses near the boundary of the cell colony, is 
given by i\ = E(l - v)h/ [Y(l + - 2v)\. 

The solution of Eq. Q with boundary conditions 
u r (r = 0) = and a rr (r = R) = can be expressed 
in terms of modified Bessel functions as 



u(r) = — (Afi 



(1 + i/)(1-2i/) 



E(l - v) 



Rh(f3r/R)A(f3), (5) 



with p = R/£ p and W)}' 1 = /3I (/3) - (^) h(P). 

As in our experiments, the resulting displacements 
and traction stresses are localized near the colony edge 
(Fig. [2|. To compare quantitatively to experiments, we 
calculate the total force, 



F(R) = 2ttF 



rdr u r (r) 



(6) 



In the large-colony limit, R ^> £ pi we find J-(R) — 
2iT(AfihR ~ R, yielding the anticipated linear growth of 



total force for large colonies. In this limit, the contractile 
active pressure dominates over internal elastic stresses 
and underlies the observed apparent surface tension. 

The theory matches the scaling of the data reason- 
ably well with £ p = 1 1 jiva and apparent surface tension 
(A/jih ~ 8 x 10 _4 N/m, as shown by the solid line in 
Fig. [3] The penetration length, £ pi is comparable to the 
spatial resolution of our measurements. For single cells, 
recent measurements have suggested apparent surface 
tensions of 2 x 10 _3 N/m in an endothelial cell [46] and 
1 x 10 -4 N/m in Dictyostelium cells [26 . From previously 
published data on a millimeter-scale adherent sheet of co- 
hesive cells, we calculated the apparent surface tension by 
integrating the average stress profiles near the sheet edge 
and found a value of about 7 x 10 -4 N/m [9]. For our cell 
colonies of thickness h ~ 0.2 /im, estimated from confocal 
imaging of phalloidin-stained colonies, the fitted value of 
the apparent surface tension implies (Afi « 4kPa. This 
value is consistent with that inferred from experiments 
in crawling keratocytes [47]. We can estimate the active 
pressure by assuming (A^i w p m /c m A m , where p m is the 
areal density of bound myosin motors, k m the stiffness 
of motor filaments, and A m their average stretch. Using 
k rn ~ lpN/nm, A m ~ 1 nm, and p rn ~ 10 3 /im -2 , we 
find CA/i^ lkPa ggHg. 

In conclusion, we demonstrate a scaling relation be- 
tween total traction force and the geometrical size of 
cohesive cell colonies adherent to soft substrates. A 
simple physical model of cohesive colonies as adherent 
contractile disks captures the essential observations and 
suggests that the apparent surface tension in the large- 
colony limit is driven by acto-myosin contractility. It is 
intriguing that a model of a cell colony with homoge- 
nous and isotropic properties is successful when the mor- 
phology of the underlying acto-myosin networks within 
the colony are patently heterogeneous and anisotropic 
(Fig. [I]). Experiments measuring the apparent surface 
tension of colonies on substrates with different stiffnesses 
and with molecular perturbations that affect the con- 
tractility of acto-myosin networks and strength of inter- 
cellular adhesions will help to illuminate the limitations 
of the current model. Additionally, the relationship be- 
tween the apparent surface tension measured here in two- 
dimensional cell colonies and the effective surface tension 
measured in three-dimensional cell aggregates [38j [40] 
needs to be established. From a cell-biology perspective, 
it will be essential to determine the molecular mecha- 
nisms that regulate a colony's apparent surface tension. 
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